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We show that the equations which describe irreversible evolution of a system can be derived from 
a variational principle. We suggest a Lagrangian, which depends on the properties of the normal and 
the so-called ’’mirror-image” system. The Lagrangian is symmetric in time and therefore compatible 
with microscopic reversibility. The evolution equations in the normal and mirror-imaged systems 
are decoupled and describe therefore independent irreversible evolution of each of the systems. 
The second law of thermodynamics follows from a symmetry of the Lagrangian. Entropy increase 
in the normal system is balanced by the entropy decrease in the mirror-image system, such that 
there exist an ’’integral of evolution” which is a constant. The derivation relies on the property of 
local equilibrium, which states that the local relations between the thermodynamic quantities in 
non-equilibrium are the same as in equilibrium. 


Most of the fundamental laws of nature are formu¬ 
lated with the help of variational principle. That is there 
exists a functional, minimization of which leads to the 
equation describing the actual evolution of the system. 
These functionals are, e.g. action in classical mechanics 
or electrodynamics, eiconal in optics, minus entropy in 
equilibrium thermodynamics. Yet, non-equilibrium ther¬ 
modynamics is probably the only area of physics, the fun¬ 
damental principle for which is formulated as inequality. 
Namely, the total entropy production, i.e. the amount of 
entropy generated in the system per unit of time, is non¬ 
negative. The entropy production can be represented as 
a bilinear combination of the flux J of material property 
and the thermodynamic force X [I|: 

a = JX> 0 (1) 

Non-negativity of the entropy production determines the 
direction of the thermodynamic process. In particular, 
this justifies the constitutive equation 

J = LxX ( 2 ) 

where lx is a so-called phenomenological coefficient. Un¬ 
like Newton’s equation of motion or Maxwell equations, 
the constitutive relation ([2]) has not been understood as 
a consequence of a variational principle. In this paper 
we formulate the variational principle, such that both 
Eq. HI and Eq. |2| follow rigorously from this procedure. 

There have been many attempts to represent irre¬ 
versible evolution of a system as being derived from a 
variational principle 3 • The difficulty lies in the fact that 
evolution is described by a parabolic differential equa¬ 
tion, which is of the first order in time. This is in contrast 
to classical mechanics or electrodynamics, the differential 
equations of which are of the second order in time. There¬ 
fore, the form of the Lagrangian for irreversible processes, 
if the one exists, is not immediately evident. If <p(r,f) is 
a variable which describes evolution of the system, e.g. 
the temperature, then the evolution equation can be writ¬ 
ten as U[(/5(r, t )] = 0, where F is an operator function. A 
common approach to construct a Lagrangian is to make it 


proportional to F. This can be either formulated explic¬ 
itly @, 3 or embedded in another framework 3, [gj. This 
approach leads to the evolution equation F[y>(r, i)] = 0 in 
a trivial way, which does not seem to add any new phys¬ 
ical understanding. It is somewhat similar to the use 
of Rayleigh dissipation function Q. Another approaches 
convert the parabolic irreversible equation to a hyper¬ 
bolic wave equation, which then implies a Lagrangian in 
a natural way. This is done by introducing a new vari¬ 
able, which satisfies a certain equati on |8l [9j] ; introducing 
non-localities in the Lagrangian [lOj, I Ilf: using restricted 
time variation 0 or fractional time derivatives 0 . 

We follow the approach of ’’mirror-image” system, 
which was originally suggested by Morse and Feshbach 
M I- This approach has not been much appreciated 
15|], probably because it was considered too artificial. 
The physical meaning of the mirror-image system, which 
gains the energy dissipated from the normal system, was 
unclear. Furthermore, Morse and Feshbach obtained the 
evolution equation only for the system with constant dif¬ 
fusion coefficient, which substantially restricted the gen¬ 
erality of their formulation. We argue that the approach 
of mirror-image system allows a clear physical interpreta¬ 
tion. In particular, it implies that increase of the entropy 
in the normal system is balanced by the corresponding 
decrease of the entropy in the mirror-image system. In 
this way we manage to derive Eq. (H| and Eq. m from 
the variational principle. In addition we derive the sec¬ 
ond law of thermodynamics for the mirror-image system. 
The important consequence of this is the understanding 
that both system should be considered in pair, while only 
one of them corresponds to a real system. Furthermore, 
we manage to derive the evolution equation and the sec¬ 
ond law of thermodynamics without restriction of the 
transport coefficients being constant, which makes the 
approach to be very general. This becomes possible due 
to employing the property of local equilibrium. 

Let the system be described by the potential ip and 
the corresponding material density p. In the case of heat 
conduction, diffusion, electric current the potential p is 
a function of (but not equal to), respectively, the tern- 
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perature T, the chemical potential pk of the component 
k, the electric potential <j>. The corresponding densities 
are the internal energy density p u , the mass density pk 
of the component k, the charge density p e . The potential 
and the density are related to each other by a material 
equation of state p = p EOS {p) . In equilibrium the poten¬ 
tial is constant throughout the system, while the density 
may vary position if there are inhomogencities. Spatial 
variation of the potential creates the flux of the corre¬ 
sponding material property, which leads to the mentioned 
irreversible phenomena. In other words, the gradient of 
the potential drives the system out of equilibrium. 

In the following analysis we shall distinguish between 
three sets of variables. The first set, p and p as de¬ 
noted above, represents the normal evolution of the non¬ 
equilibrium system. The second set, tp* and p*, repre¬ 
sents the mirror-image non-equilibrium system. It can 
be viewed as the system with ’’negative friction”, or as 
the system, where all the processes happen backwards in 
time. Finally, the third set, p eq and p eq , represents the 
equilibrium system, which does not evolve in time. Both 
p and p* depend on the position r and time t, while p eq 
is independent of either r or t. 

Following the standard variational procedure, we as¬ 
sume that there exist an action functional A[p\ = 
f Ldrdt such that the actual evolution of the system cor¬ 
responds to the extremum of this functional. We postu¬ 
late that the Lagrangian which describes the irreversible 
behavior of the system has the following form: 

L = £(p eq ) S/p -Vp* + ± ( pp * - p*p) (3) 

where £ is positive and even with respect to time reversal. 
Furthermore, £ is independent of either p or p* but may 
depend on p eq . We will discuss below the origin of such 
dependence. The important property of the Lagrangian 
(f3|) is that it is symmetric in time. If we replace t —> 
—t and simultaneously mirror-image variables with the 
normal ones, the Lagrangian will preserve its form, since 
both p and p are even in time. 

The Lagrangian © generates an infinite number of 
evolution trajectories p(r,t ) and p*(r,t). The extremal 
evolution trajectory, at which the action A reaches its ex¬ 
tremum, is the solution of the Euler-Lagrange equations 
for the Lagrangian ©: 

dL d dL dL 

-V-= 0 

dp dt dp dS7p 

( 4 ) 

dL _ d_dL_ _ dL 

dp* dt dp* dS7p* 


This is the case, in particular, for the equation of state 
PM) = P eos ’t))- 

The Euler-Lagrange equations contain derivatives of 
the Lagrangian with respect to the potentials p, p *, and 
their partial derivatives of the first order with respect 
to time, p, p*, and position, 'S/p, S/p*. The derivatives 
dL/dp and dL/dp* are evaluated straightforwardly. The 
time derivatives are evaluated in the following way. Since 
the Lagrangian © depends explicitly on the time rate of 
change of the densities, but not the potentials, we take 
the derivative of a composition of functions: 

dL dL dp* dL dp* dL dp eq 

~dip* ~ dplfy* ~ ~d/>*~dp* ~ ~dff dp eq A 

The last equality is possible due to the property of local 
equilibrium. The derivative dp eq /dp eq is independent of 
time and therefore 

d dL dp eq d dL dpi . 1 . 

dt~d/f ~ dp eq dt~d/f ~ ~dp 2 ^ ~ 2 P ^ * 

where we used the property of local equilibrium for the 
second equality. Finally, when evaluating the derivatives 
V ■ dL/dS/p and V • dL/dS/p* we use the property of 
local equilibrium one more time to write 


£(p eq ) = £{p) = £(p*) (5c) 

It is important to realize that the property of local 
equilibrium should be applied for the extremal trajectory 
onlyThis means that the above assignments should be 
made only in the evaluated expressions, i.e. after the 
variation is performed, and not in the Lagrangian © 
itself. 

As a consequence of the validity of local equilibrium 
only for the extremal trajectory, we do not have to re¬ 
quire £ to be a constant. As Lagrangian d3j suggests, 
the coefficient £ depends on the equilibrium potential p eq 
only, which is constant for the entire system. However, 
the property of local equilibrium implies that for the ex¬ 
tremal evolution trajectory £ is the same function of the 
local non-equilibrium potential p or p *. As these poten¬ 
tials may depend on position or time, so does the coeffi¬ 
cient £. 

Substituting the evaluated derivatives in the Euler- 
Lagrange equations, we obtain the following evolution 
equations: 


V • (£(p)S7p) + /> = 0 
V • {£{p*)S7p*) - p* = 0 


( 6 ) 


We assume that the extremal trajectory obeys the so- 
called property of local equilibrium. Namely, local equi¬ 
librium means that all the relations between the ther¬ 
modynamic quantities for the extremal evolution trajec¬ 
tory in non-equilibrium are valid locally in every point of 
space and time and remain the same as in equilibrium. 


The first of Eq. © represents the non-equilibrium evo¬ 
lution of the normal system, while the second of Eq. © 
represents the non-equilibrium evolution of the mirror- 
image system. 

We next consider the balance equation for the density 
in a convection-free system. It can also be viewed as 
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the balance equation in a co-moving (e.g. barycentric) 
frame of reference. The change of material property with 
respect to time is caused by the diffusional flux J. It is 
the same flux which appears in Eq. m for the entropy 
production. The balance equations for both, the normal 
system and the mirror-image system, have the same form: 


V • J + p = 0 

V • J* + p* = 0 


(7) 


The reason for this is that the diffusive flux and the time 
derivative are odd function of time, while the density is 
even. Comparing Eq. © with Eq. m we can conclude 
that 


J = £Vp 
J * = -~£Vp* 


( 8 ) 


These equations represent the constitutive relations be¬ 
tween thermodynamic forces V< p, S7p* and the corre¬ 
sponding fluxes J, J* in the normal and the mirror-image 
systems, respectively. Eq. © is equivalent to the first of 
Eq. (©. The fluxes J and J* have opposite signs. Since £ 
is positive, matter diffuses in the direction of the thermo¬ 
dynamic force in the normal system and in the opposite 
direction in the mirror-image system. 

Substituting Eq. © in the expression m for the en¬ 
tropy production (which is the same for both the normal 
and the mirror-image systems) we obtain 


binary mixture t\\ = DTpc 2 /(dpi/dci), where D is the 
diffusion coefficient and ci ,2 are the mass fractions of the 
components (ci + C 2 = 1). Finally, £ ee = a e T , where a e 
is the electric conductivity. 

An important part of the Lagrangian is the depen¬ 
dence of the coefficient t on the equilibrium potential 
p eq rather than the non-equilibrium potentials p and p*. 
As we have just concluded, the coefficient £ is the phe¬ 
nomenological transport coefficient. According to Green 
and Kubo GM2 it is equal to the integral of the cor¬ 
responding equilibrium time correlation function of the 
microscopic fluxes. It depends therefore only on equilib¬ 
rium properties of the system and should not be affected 
by the variation of the non-equilibrium potentials p and 
p* . The Lagrangian © describes therefore the linear 
reaction of the system on a thermodynamic perturba¬ 
tion. It has the same range of applicability as classical 
irreversible thermodynamics. 

According to the variational analysis [l4|, the La¬ 
grangian © implies the existence of a 4x4 tensor W 
which possesses certain symmetries. This tensor has the 
components 


W a p=L& aP -v a - - p* a — (10) 

o<PP op* p 

where p a = dp/dx a , and x a represents either one of the 
three spatial coordinates or the time coordinate, and 8 a p 
is the Kronecker symbol. The tt component of this tensor 
satisfies the following relation: 


a = £|V(^| 2 

(9) 

u* = — £ | V<^?* | 2 

We postulated that £ is positive. This means that in the 
course of evolution a may not become negative and is 
always non-negative, while cr* may not become positive 
and is always non-positive. The first of these statements 
is one of the formulations of the second law of thermo¬ 
dynamics. 

It is now evident, what is the exact form of the po¬ 
tential p for different irreversible processes. For heat 
conduction the potential is pt = 1 /T . For diffusion in 
binary mixture the potential is p\ = — (pi — P 2 )/T . For 
electric current the potential is p e = —<j>/T. Note, that 
in case of pure diffusion and electric current the tem¬ 
perature should be considered constant with respect to 
position, so that the actual driving force is —V(/zi — P 2 ) 
and — S7(j) respectively. If the temperature varies with 
position, this leads to coupled irreversible phenomena, 
which require additional treatment. 

Given this identification of the potentials, it is clear 
that the coefficient £ is the ordinary transport coefficient: 
£ qq in the case of heat conduction, l\\ in the case of 
diffusion, and £ ee in the case of electric current. These 
coefficients are related to the measured transport coeffi¬ 
cients in the standard way [l|. In particular, £ qq = kT 2 
where k is the thermal conductivity. Furthermore, in a 


iS v w “ iV=a (11) 

where the integration is performed over the volume of the 
system. We see, that W tt is the analogy of the Hamilto¬ 
nian density in classical mechanics, and Eq. m states 
that the total Hamiltonian is an ’’integral of evolution”, 
i.e. a conserved quantity. 

It follows from Eq. © that the Hamiltonian density 
corresponding to this Lagrangian is 


Wtt = £{<p eq ) v<p ■ Vp* 

(12) 

Comparing it with Eq. (© we see that 


a cr* = -W? t < 0 

(13) 


Just like Eq. it tells about the direction of non¬ 
equilibrium processes. It is evident from Eq. Cl that 
a and a* have opposite signs. This is also true for the 
total entropy productions dS/dt = f a dV and dS*/dt = 
f (7 dV. This means that entropy increase in the normal 
system is balanced by the corresponding entropy decrease 
in the mirror-image system. 

In addition, Eq. (ED gives information about the rate 
at which the entropy productions change with respect to 
time. Depending on the external conditions, the entropy 
production may have different dependence on time. In 
particular, in a stationary process in the normal system 
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the local and the total entropy productions are constant. 
In contrast, in a relaxation process in the normal system 
the entropy production decreases with time, approach¬ 
ing zero in equilibrium. It follows from the above anal¬ 
ysis that the quantity which has the density \J—a cr* is 
exactly constant irrespectively of the external conditions: 

-y- [ \]—o a* dV = 0 (14) 

at J v 

It follows, in particular, that in stationary states in the 
mirror-image system the entropy production is also con¬ 
stant (but negative). In contrast, a relaxation process in 
the normal system corresponds to a ” tightening” process 
in the mirror-image system. In this process the abso¬ 
lute value of the (negative) entropy production increases 
continuously and diverges with time. 

To understand the meaning of this relation it is use¬ 
ful to consider a simple relaxation process which occurs 
when two thermal baths which have different tempera¬ 
tures are brought in contact with each other. In this 
case the entropy is produced only at the contact of these 
baths. Volume integration in Eq. m reduces to eval¬ 
uation the value at the contact area, and instead of the 
densities a, a* one should use the total entropy produc¬ 
tions AS , AS* respectively. It follows therefore that 

AS AS* = const < 0 (15) 

Eq. (fl5l) says that the total entropy increase in the nor¬ 
mal system is inversely proportional to the total entropy 
decrease in the mirror-image system. The rate of the en¬ 
tropy change in either system is such that their product 
remains constant with time. Furthermore, the conserved 
Hamiltonian Wtt has a simple meaning for this process. 
Namely, the square of the Hamiltonian is equal to minus 
the product of the entropy productions in the normal and 
the mirror-image systems: Wj~ t = —SS*. 

The Lagrangian (j3j) itself has a simple analogy to clas¬ 
sical mechanics as well. The term —£N/p ■ 'S/ip* is analo¬ 
gous to the potential energy, while the term (1/2 ){pp* — 
p*p) is analogous to the kinetic energy. Indeed, in clas¬ 
sical mechanics if a body is not subjected to a spatially 
varying external potential, then there are no forces acting 
on it. The body remains therefore in its state, which is 
motion with a constant velocity. In non-equilibrium ther¬ 
modynamics if the potential ip is spatially constant, then 


there are no thermodynamic forces driving system out of 
equilibrium. The system remains therefore in its state, 
which is an equilibrium state. Indeed, in this case the 
evolution equations are p(r,t) = p*(r,t) = p eq = const. 
The spatial variation of the potential ip brings the system 
out of equilibrium. The term —PS/p-S/p* can be viewed 
as the potential entropy production. The analogy with 
classical mechanics is even more evident if the perturba¬ 
tion of the system from equilibrium can be described by 
an additional term to the microscopic particle Hamilto¬ 
nian, which is the case e.g. for electric current. In this 
case the potential p is proportional to the additional term 
to the microscopic Hamiltonian, which can be viewed as 
the additional energy or the entropy production which 
are potentially stored in the system. 

In conventional irreversible thermodynamics the force- 
flux relation follows from the second law of thermody¬ 
namics, which states that the entropy production is non¬ 
negative. Though the second law of thermodynamics is 
a fundamental law, it has been argued that it is not com¬ 
patible with the time reversal symmetry of the micro¬ 
scopic equations of motion. We have shown that the 
second law of thermodynamics as well as the force-flux 
relations follow explicitly from a variational procedure, 
which is compatible with the time reversal symmetry. 
Namely, that is the principle of extremal action with the 
Lagrangian given by Eq. ©■ As a consequence of the 
time reversal symmetry of the Lagrangian © , it is pos¬ 
sible to obtain two ’’second laws”: one for the normal 
system and another for the mirror-image system. En¬ 
tropy increase in the normal system is balanced by the 
entropy decrease in the mirror-image system. 

An important consequence of the suggested variational 
approach is that the evolution equations © for normal 
and mirror-image systems are decoupled. This means 
that local evolution of the normal system is independent 
of the local evolution of the mirror-image system. One 
can say that the two systems evolve in separated time¬ 
lines without affecting each other. The sign of the en¬ 
tropy production is determined by the timeline to which 
the system belongs. This sign cannot change during evo¬ 
lution. It is always positive for the normal system and al¬ 
ways negative for the mirror-image system. One can also 
formulate the inverse statement: the direction of time is 
determined by the sign of the entropy production. 
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